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The same light-front wave functions of the proton are involved in both the anomalous magnetic 
moment of the nucleon and the Sivers function. Using the diquark model, we derive a simple 
relation between the anomalous magnetic moment and the Sivers function, which should hold in 
general with good approximation. This relation can be used to provide constraints on the Sivers 
single spin asymmetries from the data on anomalous magnetic moments. Moreover, the relation 
can be viewed as a direct connection between the quark orbital angular momentum and the Sivers 
function. 
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The quark orbital angular momentum [ij (or quark transverse motion) plays an important role for understanding 
the spin and quark structure of the nucleon, since as shown by many studies 0, H, 0, [H> H 0] , it is the missing block of 
the total nucleon spin. Also many interesting phenomena or observables require the presence of quark orbital motion, 
among which the Sivers function Q has attracted a lot of interest, since it is an essential piece in our understanding 
of the single spin asymmetries (SSA) observed in semi- inclusive deeply inelastic scattering (SIDIS). These SSAs have 
been measured recently by both the HERMES @,H3| and COMPASS [HHl] Collaborations. Denoted as f^ T (x,k\), 
the Sivers function describes the unpolarized distribution of the quark inside a transversely polarized nucleon, which 
comes from a correlation of the nucleon transverse spin and the quark transverse momentum. Although this is a T-odd 
type correlation, it has been found that final state interaction [l3| (FSI) between the struck quark and the spectator 
system can produce the necessary phase for a non-zero Sivers function, and its QCD definition [3, EH [H, [TtJ has 
just been worked out. Besides the single spin asymmetry, another important feature of the Sivers function is that it 
encodes the parton's orbital angular momentum (L z ) inside the nucleon. This comes from the fact that the Sivers 
function requires the nucleon helicity to be flipped from the initial to the final state, while the quark helicity remains 
unchanged. A convenient tool to study this kind of single spin asymmetry (or the corresponding Sivers function) is the 
light-front formalism [18 1, which can express the Sivers function as the overlap integration of light-front wavefunctions 
differing by AL Z = ±1 pj|[3. 

The same kind of overlap integration [20L l2ll H2] of light-front wavefunctions (with 
J z = ±1/2 in the initial and final states) also appears in the anomalous magnetic moment of the nucleon, which 
apparently encodes the quark orbital angular momentum [211 ]. Therefore, it is interesting to find relations between 
the Sivers function and the anomalous magnetic moment of the nucleon, which is the main goal of this work. With 
such a relation one can constrain the Sivers function and the related asymmetries from data on nucleon anomalous 
magnetic moments, and viceversa. Also, the relation can be viewed as a direct connection between the quark orbital 
angular momentum distribution and the Sivers function. 

The proton state can be expanded in a series of Fock states \n, x^P + , kx, Xi) with coefficients ip n / p (xi, k±, Ai), which 
are the light-front wavefunctions of the proton: 

* p (P+,P-,Ox) = Y,M^x,Xi)\n,XiP+, k ± ,\i). (1) 

n 

k+ . . 

Here Xi = -pt? is the light-front momentum fraction of the quark, Ai denotes the helicity, and its transverse 
momentum. The wavefunctions are Lorentz-invariant functions of the kinematics of the constituents of nucleon, X{ 
and k_u, with YTi x i = 1 and Yh k -L* = °-L- 

As pointed out before, the Sivers function requires that the nucleon wavefunctions in the initial and final state 
differ by one unit of orbital angular momentum, and final state interactions play a crucial role. It describes the 
interference of two amplitudes which have different initial proton spin J z — ±i but couple to the same final state: 
lm[M*(J z — +1/2 — > F)M(J Z — —1/2 — > F)]. This can be realized by a gluon exchange between the struck quark 
and the spectator system. There have been already attempts 0, using the proton light-front wavefunctions, to 
find a formula to calculate the Sivers functions. In those papers the final state interaction phase needed for Sivers 
functions has been incorporated into the wavefunctions. Another possibility is to express the Sivers function as the 
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product of wavefunctions and the final state interactions term: 

ft T ex ^^(^^x^^G^.^.k^.kl^^.k'x.AO, (2) 

n 

where G(x' i ,'k' ± i ,Xi,k±,i) is the final state interaction term, and the light-front wavefunctions in this equation are 
the usual ones which do not contain the final state interactions phase. 

In the scalar diquark model, the calculation gets simplified and some instructive result can be derived. In this 
model the kinematics is determined by the momentum fraction and the transverse momentum of the struck quark 
(x, kjj, since those of the spectator diquark can be related to (x, kj_) by xd — 1 — x, k±£> = kj_. The advantage of 
the diquark model is that it is simple, and it works quite well in describing other properties of the proton, such as the 
helicity distribution [24[, the transversity [25| and the form factors [26| of the proton. The light-front wavefunctions 
of this two-body Fock state have the simple forms [13, [2fJ : 

/ T 2 I 1 \ ( + fe 1 +ife 2 ) / i \ (3) 



for J z = +1/2 and 



^(x,k ± ) = <±^%(x,k ± ) , 



for J z = -1/2. Here 



ufxk,)- g/vT 3 ^ , 5) 

^> K Xj M 2 -(k x +m 2 )/x-(k ± +M 2 )/(l-xy { > 



and M, m and Md are the masses of the proton, the quark and the diquark, respectively. 

According to Fig. [U a formula to calculate the Sivers function can be given by the overlap integration of the proton 
light-front wavefunction: 



uL p fly i i 

^f#(x,k x ) = l J2J ^^e(^,k i ,A l )G(x,x',k i ,k 1 )^(x',k' ± ,AD, (6) 



where k L = k 1 — ik 2 , and G is the kernel which contain the final state interaction. In the scalar diquark model, the 
kernel has a simple form [l3|, H3| : 

G(x x 1 k , kM - iCpa ^ x ~ x ') (71 

which can be calculated in the eikonal approximation. Final state interactions arise from gluonic exchanges between 
the struck quark and the spectator system, which are necessary in order to insure gauge invariance of the kj_ -dependent 
distributions. In our case we use a o ne-g luon exchan ge ( with momentum k — k') approximation, which is commonly 
used in most model calculations [H, [TH, [H, [2iJ l30ll3ll |32| . Here A g is the mass of the gluon, which is needed in 
general to regularize the infrared divergence in the integration, while in our case it can be set to A g — > 0. The sign 
in Eq. (J7J) is according to the Trento conventions [33| . With Eqs. ([3]), (0J and ©, one can obtain the result for 
Sivers function that has given in An expression similar to ([6|) has been given in Ref. where the Sivers 

function (or the SSA) is expressed as a product of a generalized parton distribution (GPD) and a FSI term in impact 
parameter space, while our result is expressed in momentum space, which can be connected with the result in Ref. [27| | 
by Fourier transformation. 

The same set of wavefunctions (spin-flipped in the initial and final states) appearing in @ also appears in the 
expression of the Pauli form factor or anomalous magnetic momentum of the proton. In light-front formalism the 
Pauli form factor is identified from the helicity-flip vector current matrix elements of the J + current plj 



J+(0) 



2P+ 



3 




FIG. 1: Diagram to calculate the Sivers function. The arrows T / I show the polarizations of the nucleon. 



Furthermore, from Fig. [2| one can express the Pauli form factor in terms of the light-front wavefunctions as 

F2{q 2 ) _ \- f d 2 k ± dx 



n j 



2M 
where 

= k ± + (1 - x)q ± , (10) 

here g is momentum of the virtual photon, the momentum transfer between the initial and final nucleoli. 

The anomalous magnetic moment is defined from the Pauli form factor in the limit of q = 0: k = (e/2Af)i 7 2(0). 
After some algebra calculation, the anomalous magnetic moment can be expressed in terms of a local matrix element 
at zero momentum transfer: 

j n ijtj ^ ' 

From ^ and (fTTj) we see that the same set of light-front amplitudes, with the orbital angular momenta differing by 
AL Z — ±1 between the initial and final state, appears in the calculation of the Sivers function and the anomalous 
magnetic moment. Recently a study [35j has shown how the non-vanishing anomalous magnetic moment constrains 
the quark orbital angular momentum. 

It is clear from the previous expressions that the relation we are seeking between the Sivers function and the 
anomalous magnetic moment will necessarily have to be approximate. This can be also seen from the fact that the 
Sivers function, although leading twist, has (Q 2 ) scale evolution dependence, which is not the case for the anomalous 
magnetic moment. 

In the scalar diquark model, the calculation of the anomalous magnetic moment at the quark level has been given 
in [2l| as 

k \ -> f d k±dx T , , , , . , / d id \ , , , , . , . 

n=±l/2 J V 7 

Using the wavefunctions given in ([3]) and (j4|), one gets the result 

Mg 2 f l (l-x) 2 (m + xM) 



(1 — x)m 2 + xm d — x{l — x)M 



where we have defined: 



A = ^(l-x)(xM + m); ( 1.4 ! 



r 2 

B = (1 - x)m 2 + xM 2 d - x(l ~ x)M 2 . (15) 
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FIG. 2: Diagram to calculate the Pauli form factor (or anomalous magnetic moment). The arrows j / J, show the polarizations 
of the nucleon. 



One can calculate the lowest kj_-moment of /{ip 9 (x, kj_) defined as (A 9 = 0) 

f^(x)=Jd 2 k x f^(x,k 2 x ). (16) 
According to Eq. |6j), which gives /{^ 9 (a;, kj_), and the wavefunctions given in (J3j> and {4|, we directly get 

ftA*) = (17) 



Defining 



which satisfies the normalization condition 



(1 - x)A 

Kq{x) = B ( 8) 



l 

dxK q (x) = Kg, (19) 



we arrive at a simply relation between n q and f^- q (x) in the quark-diquark model 



/rr!^-^^^,-). ,20) 



The term 1/(1 — x) in this equation is consistent with the result in Ref. [19}. From this we can get: 

/ dx(l - x)f^ T q (x) = - [ dx - C ^ as n q (x) 
Jo Jo iz 

n 2 CFa s 

12 Kq ' 

An approximate relation can also be obtained from the above equation, as 



(21) 



(1 - x)Ffr = ~^^Kg, (22) 

where F^ q is the first x- moment of f^ q (x) denned as F^j, q = f Q dxf^ q (x) , and (1 — x) is the average value of 
(1 — x) in the proton. Therefore we obtain a interesting relation that provides a constraint on the Sivers function 
from anomalous magnetic moment data. It suggests that the moment of the Sivers function is proportional to the 
anomalous magnetic moment contributed by the same quark. Although this relation as been obtained in the diquark 
model, it should hold in general with good approximation. This can be seen from the general relations for the 
anomalous magnetic moment and for the Sivers function in terms of light-front wavefunctions, as given by Eqs. ^ 
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and (jlip . respectively. Since both quantities calculated from the same wavefunctions, and these are the ones that 
vary more strongly in the range of values covered by the integrals, it is appropriate to apply the mean value theorem, 
and therefore get the same type of relation as the one we have found here. In this sense our approximate result can 
be considered more general than the diquark model, and thus be applied to each quark individually. 

A more fundamental object than the anomalous magnetic moment is the spin-flipped generalized parton distribution 
(GPD) E(x,£,t) @,[3i,[33,[3i|. Actually the function n q (x) is E(x,£,t) at the forward limit (£ = and t = 0): 

K q (x) = E q (x,0,0)- (23) 

Thus in the scalar diquark model, f^, q (x) is proportional to E q (x, 0, 0). According to Ji's sum rule [4( ( J q is the total 
angular momentum carried by quark flavor q): 

J dxx(H q (x,Z,t) + E q (x,£,t)) = \JS)> (24) 

which also holds in the forward limit. We see that the Sivers function is related to the angular momentum of the 
parton inside the nucleon, and therefore it is in fact sensitive to the orbital angular momentum of the quark. There 
is then the possibility to get information of the quark orbital angular momentum from the Sivers functions. 

Although the relation given in (j2"0|) and (|21[) is a simple result based on the approximation of the scalar diquark 
model, we can still apply the relation to given some prediction on the Sivers asymmetry of the meson production in 
SIDIS processes, such as the ratio of the asymmetries between different final mesons, since in this case the model 
dependence is reduced. The Sivers asymmetry is proportional to 



ASiv (So e cJlT D ll / 9t -\ 



which can be extracted by the weighting function sin(</> — 0s), here <j) and (f>s denote respectively the azimuthal 
angles of the produced hadron and of the nucleon spin polarization, with respect to the lepton scattering plane, D\ 
is the unpolarized fragmentation function. We will focus on the large z regime that the valence quark contribution 
dominates, and the disfavored fragmentation function can be ignored. 

Since the Sivers function and the anomalous magnetic moment "share" the same set of the proton wavefunctions, 
as shown in Figs. [T] and [21 one can start from the data of the anomalous magnetic moment to provide constraints on 
the proton wavefunctions, and then on the Sivers function. Similar methods have been used in Ref. [l9j , where a small 
Sivers asymmetry on a deuteron target has been predicted, and Ref. [39| . where the sign of the Sivers asymmetries 
for different hadron targets combining different fragmenting hadrons has been classified. 

As figured out in Ref. [l9l |. the quark contribution dominates over the gluon contribution in the case of Sivers 



functions [40j, which is also the result of an argument based on large N c considerations [4l|. There are also phe- 
nomenological supports of this conclusion from the SIDIS experiment from COMPASS of CERN [ll[, as pointed out 
in Ref. [l3 | , and the analysis on hadron production of n given in Ref. [42| . Therefore in this work we only consider 
the quark contribution to the Sivers functions and the corresponding asymmetry. 

One can constrain the proton wavefunctions by normalizing each u and d quarks contributions to the anomalous 
moments n p — 1.79, K n — —1.91. Isospin symmetry implies 

K d/n — K u/p, (26) 
«u/n = K d/p- (27) 

In the valence quark approximation we have: 

k p = (2)(2/3)k„ /p + (-l/3)K d/p , (28) 
n n = (2)(-l/3K /p +(2/3)/« d/p . (29) 

Thus one has n u / p — 0.835, n d / p = —2.03. In the following we use n u and Kd to represent n u / p and Kd/p, respectively. 
Then we can write the ratio of the asymmetries between 7r + and 7r~ at large z: 



= -3.3. (30) 
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Also we have 



2e d n d 



e * Kd = -1.15, (31) 



6.6. (32) 



For the above result we used isospin symmetry for the quark fragmentation functions: 

Df /u = Df /d = 2D? /v = 2D? /d , (33) 
Df /U = 2Df' d . (34) 

The results show that in the large z region, the Sivers asymmetry of 7r + is 3 times larger than that of 7r - and with 
opposite sign, which is consistent with the recent HERMES result [l0( where at z ~ 0.6 a four times larger asymmetry 
of 7r + is measured; the asymmetries of ir° and ir~ are similar in size; the asymmetry of K + is much larger than that 
of K , nearly one order of mag nitude; and Af%(K°) ~ since k s = in valence approximation. 

In summary, both the formula for calculating Sivers function and that for calculating the anomalous magnetic 
moment of the proton, can be expressed in terms of the same set of the light-front wavefunctions, with hclicity flipped 
from initial state to final states. Using the overlap representations of both Sivers functions as well as the anomalous 
magnetic moment, we give a simple relation between these two observables, in the approximation of the scalar diquark 
model. This relation is applied to provide constraints on the Sivers single spin asymmetries in the valence regime from 
data on anomalous magnetic moments. Also, the relation can be viewed as a connection between the quark orbital 
angular momentum and the Sivers function. 
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